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ABSTRACT 

We report on the results of a Submillimeter Array interferometric observation 
of the proto-planetary nebula CRL 618 in the ^^CO J = 6—5 line. With the new 
capability of SMA enabling us to use two receivers at a time, we also observed si- 
multaneously in the ^^CO J = 2-1 and ^^CO J = 2-1 lines. The ^^CO J = 6-5 
and ^'^CO J = 2—1 lines were first interferometrically observed toward CRL 618. 
The flux of the high velocity component of the ^^CO J = 6—5 line is almost fully 
recovered, while roughly 80% of the flux of the low velocity component is resolved 
out. The low recovery rate suggests that the emission region of the low velocity 
component of the ^^CO J = 6—5 line is largely extended. Continuum emission is 
detected both at 230 and 690 GHz. The flux of the 690 GHz continuum emission 
seems to be partially resolved out, suggesting dust emission partly contaminates 
the 690 GHz continuum flux. The cavity structure, which has been confirmed in 
a previous observation in the ^^CO J = 2—1 line, is not clearly detected in the 
^^CO J = 6—5 line, and only the south wall of the cavity is detected. This result 
suggests that the physical condition of the molecular envelope of CRL 618 is not 
exactly axial symmetric. 
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Introduction 



Stellar evolution from the asymptotic giant branch (AGB) phase to the planetary neb- 
ula (PN) phase is very rapid. This transient phase between the AGB and PN phases is 



often called the pr oto-planetary nebula (PPN) phase and/or post- AGB phase ( lKwoklll993 



van Winckelll2003l ). PPNe are considered to play an important role in a wide variety of as- 
trophysical problems: for example, the shape and shapi ng of PNe, chemical evolution of th e 
evolved stars, and synthesis of organic matter in space (IBalick fc Frankll2002l : lKwokll2004j ). 
Unfortunately, however, the number of PPNe identified is rather limited mainly due to their 



transient time (see, e.g.. ISzczerba et al.ll2007l ). 



In the last few decades, the well-known example of a PPN, CRL 618 (= RAFGL 618 
= IRAS 04395+3601 = Westbrook Nebula), has provided a unique opportunity to investigate 
the nature of at least one particular case of PPNe. CRL 618 entered its PPN phase about 
100 years ago (IKwok fc Bignelllll984l ). and the central BO star is surrounded by a compact 
Hn region (its angular size is roughly 0.2^^x0.4^0 visible through cent i meter- and millimeter- 
wave continuum emission (IWvnn- Williams! 119771 : iKwok fc FeldmanI Il98ll : iKwok fc Bignell 
1984J : iMartin-Pintadd Il993l ) . The flux of the free- free emis sion at A = 1.3 and 3 mm has 



been continuously changing over the last 25 years (see, e.g.. ISanchez Contreras et al.ll2004j ). 
implying the rapid evo lution of this object. The envelope of CRL 618 includes a rich set of 



molecular species (e.g., iBuiarrabal el at 



Cernicharo et al. 2001 



200^ 



1988 



Cernicharo et al. 1989: Fukusaku et al. 1994 



2; Remijan et al. 2005[). and has been repeatedly observed in 



molecular radio lines, especially that of CO (e.g..lShibata et al.lll993l : lYamamura et al.lll994 
Hajian et al.lll996l : ISanchez Contreras et al.l 120041 ). 



Through radio observations in molecular rotational lines, four different kinematical com- 
ponents have been identified in the molecular envelope of CRL 618. In Figure 1, we present 
a schematic view of the molecular envelope of CRL 618. The first component (corresponding 
to (1) in Figure 1) exhibits an almost round shape in CO maps, and its angular size is larger 
than 20". The expanding velocity of this spherical component is roughly 17.5 km s~^, and is 
interpreted as the remnant of the spherical mass-loss process of the central star during the 
AGB phase. The second component (corresponding to (2) in Figure 1) exhibits a bipolar 
flow observed up to axial distances of ±2.5" from the nebula center. The expansion velocity 
of this component reaches up to 350 km s~^. Incidentally, the prope r moti on of this high 
velocity bipolar flow has been measured by ISanchez Contreras fc Sahail (120041 ). and the value 
0.045" yr~^ corresponding to a distance of 900 pc if we assume the inclination angle 



IS 



of the p olar axis is 30°: this distance is consistent with that estimated by an independent 
method (lGoodrichlll991l) . An a dditional two kinematical components were recently found by 
Sanchez Contreras et al.l (120041 ) in their high-resolution (~ 1") interferometric radio obser- 
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vation: an extended structure elongated in the polar direction up to an axial distance of ±6" 
from the nebula center, which is expanding at a velocity of about 22 km s~^ (corresponding 
to (3) in Figure 1) and a dense, inner torus-like core expanding at a velocity less than 12 
km s~^ (corresponding to (4) in Figure 1). 

Although numerous observational and theoretical efforts have been made on CRL 618, 



tectable in the HCO+ J 



tin the CO J = 2-1 line ( 


Sanchez Contreras et al. 


Sanchez Contreras et al. 


(2004 


) suggested that a 



passage of a shock across the torus-like structure plays an important role in explaining 
this phenomenon. Second, the light curve of t he intensity of the free-free em ission ex- 
hibits a sudden increase over the past 5 years (jSanchez Contreras et al. 2004), implying 



that CRL 618 exhibits the a c tivity of a pqst-AGB wind, which is ve r y poorly understood 



(e.g., iBujarrabal et al.l 1 19941 : lJura fc ChenI l2000l : ISteffen et al.l 1200 ll : ICarcia-Segura et al. 
20051 ). Third, the molecular envelope of CRL 618 exhibits a hyb rid chemistry simulta- 



neously including both oxygen- and carbon-rich environments (e.g., IBujarrabal el at.lll988 



Cernicharo et al.lll989l ). Such a mixture of the C/O-rich chemistry has been r ecognized in 



some late-type stars (e.g . CRL 2688, OH 231.8+4.2 and IRAS 19312+1950, 



Lucas et al. 



19861 : Morris et al.l 119871 : iNakashima fc Deguchil l2000l . I2004J . l2005l : iMurakawa et all 120071 ). 
but the origin of the mixture is still unclear. 

Further observations of CRL 618, especially using new instruments with new capabilities, 
would be valuable in understanding the nature of this interesting object. In this paper, 
we report the first results of a Submillimeter Arra}0 (SMA) interferometric observation of 
CRL 618 in the ^^CO J = 6—5 line using the new capability of SMA enabling us to observe 
at 690 GHz. This observation has been made as a part of the SMA 690 GHz observation 
campaign. In addition to the ^^CO J = 6—5 line, we briefly report the results of observations 
in other molecular lines including the ^^CO J = 2—1 and ^^CO J = 2—1 lines. The outline of 
the paper is as follows: in Sect 2, details of the observations and data reduction are presented. 
In Sect 3, results of the observation are presented, including spectra and intensity maps of the 
CO lines, intensity and position measurements of the continuum emission and identification 
of the detected molecular lines. In Sect 4, we discuss the properties of the ^^CO J = 6—5 
line. Finally the results of the present observation are summarized in Sect 5. 



^The Submillimeter Array is a joint project between the Smithsonian Astrophysical Observatory and the 
Academia Sinica Institute of Astronomy and Astrophysics, and is funded by the Smithsonian Institution and 
the Academia Sinica. 
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2. Details of Observations and Data Reduction 

Interferometric observations of CRL 618 were made with SMA on February 18, 2005. 
Data were obtained under very good atmospheric conditions with a zenith optical depth at 
230 GHz of 0. 05 o n average. The instrument has been described in detail by lHo et al.l (120041 ) 



and lOhashil (120051 ). With the capability of SMA enabling us to use two different receivers at 
a time, we observed simultaneously in four different frequency bands: 690.572-692.549 GHz, 
680.572-682.550 GHz, 219.598-221.575 GHz and 229.600-231.576 GHz. These frequency 
ranges cover three CO lines: the ^^CO J = 6-5 , ^^CO J = 2-1 and ^^CO J = 2-1 lines. 
The array used 6 elements (with a diameter of 6m) in the compact array configuration. The 
baseline length ranged from 10 to 70m. The field of view of a single antenna was 18" and 54" 
at 690 and 230 GHz, respectively. The observation was interleaved every 15 minutes with 
nearby gain calibrators, 3clll and Titan, to track the phase variations over time. Because 
the angular size of Titan was roughly 0.88" x 0.68", Titan was not an ideal point source 
compared with the small size of the beam (~ l"-3") of the present observations. Therefore, 
in the gain calibration, we used the "non-point" flag of the "gain-cal" command in the IDL- 
MIR package, which enabled us to use a planet /satellite as a calibrator (IDL-MIR is a data 
reduction software developed by the SMA project). The non-point flag function assumes the 
calibrator is a uniform disk; the size and flux of the calibrator are automatically retrieved 
and calculated by the internal theoretical calculation. The absolute flux calibration was 
determined from observations of Ganymede, and was approximately accurate to within 15%. 
We adopted the antenna-based solution in the pass band calibration of the 690 GHz data. 
With enough signal to noise ratio on each channel, the antenna-based bandpass solution is 
robust and consistent with baseline-based calibration results. For the 230 GHz data, we 
adopted a standard calibration method using the baseline-based solution. The final map 
has an accumulated on-source observing time of about 4 hours. The single-sideband system 
temperature ranged from 2000 to 3000 K at 690 GHz and from 150 to 200 K at 230 GHz, 
depending on the atmospheric conditions and the telescope direction. The SMA correlator 
had a bandwidth of 2 GHz with a resolution of 0.812 MHz. The velocity resolution was 
0.35 km s~^ at 690 GHz and 1.0 km s~^ at 230 GHz. The phase center of the map was 
taken at R.A.= 04M2"53.67^ decl.= 36°06'53.2" (J2 000). Image pro ce ssesing of the data 



was performed with the MIRIAD software package (ISault et al.l Il995l ) . We used robust 
weighting for the 690 GHz data. The robust weighting of the visibility data, which is an 
optimized compromise between natural and uniform weighting, gave a 1.0" x 0.8" CLEAN 
beam at 690 GHz with a position angle of —12°. For the ^^CO and ^^CO J = 2—1 lines, a 
super-uniform weighting of visibihty data (with a suppression region of 8" x 8") was adopted 
to enhance the angular resolution, and this weighting gave a 2.6" x 2.4" CLEAN beam with 
a position angle of +11.2°. The width of the bands for the continuum maps at 690 and 230 
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GHz were roughly 1.5 and 1.9 GHz, respectively. 



Results 



3.1. Spectra of CO lines 



Figure 2 shows spatially integrated spectra of the observed CO lines. The line emission 
was integrated over a circular region with a diameter of 15". The continuum flux was calcu- 
lated in the line-free channels, and the continuum emission was subtracted from the spectra 
to retrieve only the line emission. As reported in previous papers, the ^^CO J = 2—1 
line clearly shows the two different compor ients (i.e., high veloci t y pedestal component 



and low velocity strong component; see , e.g.. iBuiarrabal el at 



Cernicharo et al. 19891: Neri et al. 1992 : Meixner et al. 



1998 



1988 



Buiarrabal et al 



Sanchez Contreras et al. 



2001 



2004 



Sanchez Contreras fc Sahaill2004j ) . In comparison with a single dish spectrum (IBujarrabal el at 
19881 ). there is no significant loss of the flux in the high velocity component of the ^^CO 
J = 2—1 line, but roughly 45% of the flux of the low velocity component is resolved out. 
In fact, the line profile of the high velocity component of the ^^CO J = 2— 1 line is quite 
similar with that obtained in previous single dish observations, but the peak intensity of the 
low velocity component is clearly lower than values of the single dish measurements. 

The ^^CO J = 6—5 and ^^CO J = 2—1 lines have been interferometrically observed 
for the first time, and, at a glance, these lines exhibit different line profiles compared with 
the ^^CO J = 2—1 line. The high velocity component of the ^^CO J = 6—5 line is clearly 
detected, exhibiting the maximum line width of about 350 km s^^ at the intensity level, 
while the low velocity component is relatively weak. This weak intensity seems to be caused 
by missing flux resolved out by the interferometry. In fact, t he peak intensity o f the ^^CO 
J =6-5 line measured by a single dish observation is 320 Jy (IHerpin et al.ll2002l ). while the 
peak intensity of the ^^CO J =6-5 line measured in the present observation is 55.3 Jy. In 
comparison with the Herpin's single dish observation, roughly 80% of the flux of the low 
velocity component of the ^^CO J =6-5 line is resolved out in the present observation, 
even though there is no significant missing flux in the high velocity component of the ^^CO 
J = 6—5 line. The large amount of missing flux of the low velocity component implies 
that the emission region of the low velocity component of the ^^CO J = 6—5 line is largely 
extended. 

The ^■^CO J = 2—1 line exhibits only the low velocity co mponent, presumably r elated to 
the sensitivity. In comparison with a single dish observation (IBujarrabal el at.lll988l ). there is 
no significant flux loss in the low velocity component of the ^^CO J = 2—1 line. In all three 
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spectra of the CO lines, we can see an absorption feature at Visr ~ —40 km s~^. Although it 
is somewhat difficult to see the absorption feature in the ^^CO J = 2—1 line (at least, there 
is no absorption below the continuum level), the line profile of the ^^CO J = 2—1 line is 
certainly asymmetric: it shows a weaker blue side compared with the red one (this is clearly 
seen in the p—v diagram given in a later section; see. Figure 8). Therefore there is in fact 
absorption of the blue side of the low velocity component in the ^^CO J = 2—1 line. The 
profile of the low velocity component showing a spiky feature is somewhat different from that 
of the ^^CO J = 2— 1 line showing a parabolic profile. This difference is presumably due to a 
relatively thin optical depth of the ^^CO J = 2—1 line. The spiky profile possibly suggests a 
complex kinematics in the inner part the molecular envelope. In the bottom panel of Figure 
2, we can see a broad emission feature at Visr ~ —300 km s~^, which can be explaned by 
leakage of the strong ^^CO J = 2—1 line lying in the other side band. This conclusion is 
based on the fact that he peak velocity of this feature exactly corresponds to the backside 
of the ^^CO J = 2—1 line, the line profile is very similar with that of the ^^CO J = 2—1 
line, and also because no possible lines are found in molecular line catalogs. (Note that this 
leakage problem has already been fixed, and the current SMA system has no problems even 
in the case of very strong lines such as CO lines.) 



3.2. Intensity maps of CO lines 

Figure 3 shows the velocity channel maps of the ^^CO J = 6—5 line. The background 
gray scale represents the continuum emission at 690 GHz, and the origin of the map is taken at 
the intensity peak of the continuum emission. In Figure 3, we can cl early see a velocity gradi 
ent in the E-W direction as well as the ^^CO J = 2—1 map given by lSanchez Contreras et al. 



(2004). Interestingly, we cannot see the cavity structure found by ISanchez Contreras et al. 
( 120041 ) in the ^^CO J = 2—1 line, although the angular resolution of the present observation 
in the ^^CO J = 6—5 line roughly equals that of Sanchez Contreras's observation in the 
^^CO J = 2—1 line. Both ends of the high velocity component (i.e., from —166 to —116 
km s~^ and from 35 to 115 km s~^) are clearly extended in the E-W direction up to roughly 
±2.5" from the nebula center. This size of the extension in the E-W direction is the same 
as that seen in the ^^CO J = 2—1 line. 

In Figure 4, we present the velocity integrated intensity maps of the ^^CO J = 6— 5 



line superimposed on the HST WFPC2 Ho+continuum image (ITrammell fc Goodrichll2002l ). 
The cross indicates the intensity peak of the radio continuum emission at 690 GHz. The 
direction of the elongation of the ^^CO J = 6—5 feature coincides with that of the optical 
feature seen in the HST image. The emission region of the high velocity component of the 
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^^CO J = 6—5 line is limited to the op tically faint region (i.e. the vic inity of the nebula 
center). In comparison with the maps in lSanchez Contreras et al.l (120041 ) . the spatial size of 



the emission region of the low velocity component of the ^^CO J = 6—5 line is smaller than 
that of the ^^CO J = 2— 1 line. This small size is due presumably to the large missing flux. 

Figure 5 shows the channel velocity maps of the ^^CO J = 2—1 line. In contrast to the 
^^CO J = 6—5 line, the velocity gradient is not clear. This is due to a low signal-to-noise 
ratio in the high velocity component of the ^^CO J = 2—1 line. In Figures 6 and 7, we present 
the velocity channel maps of the ^^CO J = 2—1 line. We present the maps of the low velocity 
component in Figure 6, and also present those of the high velocity component in Figure 7. 
The features of the ^^CO J = 2— 1 line seen i n Fig ures 6 and 7 are in good agreement with 
the maps given by ISanchez Contreras et al.l (120041 ): we can clearly see a velocity gradient 
in the E-W direction. A cavity structure seen in Sanchez Contreras's maps is not clear in 
Figure 6. This is due to the low angular resolution of the present observation. In Figure 7, 
both ends of the high velocity component are extended to in the E -W direction up to ±2.5^ ^ 



from the nebular center. This is also consistent with the results of Sanchez Contreras et al. 



( I2OO4J ) . We cannot see the spherical halo seen in t he Sanchez Contreras's naaps du e to a lack 
of sensitivity. In fact, a 1 a level in Figure 2b in ISanchez Contreras et al.l (120041 ) is 45 mJy 
beam~^ (beam size is 1.1" xO. 9"); in contrast, a 1 ex level in Figure 6 is 100 mJy beam~^ 
(beam size is 3.3" x 2.8"). Missing flux, in principle, could be a reason for the absence, but 
we recovered 55% of the flux of the low velocity component, and this recovery rate is larger 
than the 40% of Sanchez Contreras's observation. 

For a better understanding of kinematical properties, we present position- velocity (p—v) 
diagrams in Figure 8. Th e present data are superirnposed on the p—v diagram of the ^^CO 
J = 2—1 line taken from ISanchez Contreras et al.l (120041 ). The left and right columns rep- 
resent p—v diagrams in the axial and perpendicular directions, respectively. In the middle 
two panels, we can see that the emission of the ^^CO J = 2—1 line is not detected in the 
high velocity ends due to a lack of sensitivity. The feature of the ^^CO J = 6—5 line is 
basically consistent with that of the ^^CO J = 2—1 line, even though the spatial size of the 
low velocity component of the ^^CO J = 6—5 line is smaller than that of the ^^CO J = 2—1 
line. 



3.3. Continuum emission 

The continuum emission of CRL 618 was detected both at 230 and 690 GHz. The 
continuum emission at 690 GHz has been interferometrically observed for the first time. 
These continuum emissions seem to exhibit a point-like feature in the intensity maps, but 
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the Mf-distance versus amplitude plot of the 690 GHz continuum emission shows that the 
amplitude increases with decreasing Mt>-distance (although we do not present the plot here), 
suggesting the 690 GHz continuum flux is partially resolved out. The integrated fluxes of 
the continuum emission are 2.0 and 2.7 Jy at 230 and 690 GH z, respectively. The 230 



GHz continuum flux is consistent with the latest measurement by ISanchez Contreras et al. 



( I2OO4J ). The spectral index calculated by the 230 and 690 GHz measurements is i/"'^. This 
is not consistent with the spectral index expected from optically thin free-free emission 
(z/^°-^). Inci dentally, the turnin g-off point of the spectrum of the central HII region is roughly 
at 45 GHz (IKnapp et al.lll993l ). Because the angular size of the central ionized region is 
0.2" xO. 4", the free-free emission is not expected to be resolved by our synthesized beam. 
Therefore, the extended feature in the 690 GHz continuum emission presumably originates 
in the dust component in the envelope. The positions of the intensity peak of the continuum 
emission measured by two-dimensional Gaussian function fltting are 4'^42™53.58'^, 36°06'53.4" 
(J2000) at 230 GHz (FWHM and position angle are 3.27"x2.84" and -19.2°, respectively), 
and 4M2"^53.59^ 36°06'53.3" (J2000) at 690 GHz (FWHM and position angle are 1.07"x0.88" 
and of —13.0°, respect ively). These position s coincide well with that determined by a low 
frequency observation (jWynn-Williamslll977l ) . but are slightly shifted from the intensity peak 
of the CO J = 6—5 emission (see lower panel in Figure 4. We will discuss this matter later 
in Sect 4). 



3.4. Other detected lines 



In Figure 9 we present a spectrum including all lines detected in the 230 GHz bands. The 
continuum emission was subtracted in Figure 9. In total we have detected 36 emission lines 
above the 5 a level, and 24 out of the 36 detected lines are identifled to known molecular lines. 
We used 'WinSpectra" , the line identiflcation software developed by the ODIN group at the 
University of Calgary to i dentify the lines. W inSpectra uses the JPL and CDMS catalogs 
(IMuller et al.l I2OOII . l2005l : iPearson et al.l |2005| ) in a built-in form to retrieve the molecular 
data. The results of the line identiflcation are summarized in Table 1. Most lines seen 



in Figure 9 exhibit a P Cygni proflle as reported in previous observation s (jBujarrabal el at. 



19881 : IWvrowski et al.ll2003l : IPardo et al.ll2004l : IPardo fc Cernicharoll2007l ) . This suggests that 
most of t he detected lines origin ate from an expanding and accelerating molecular envelope 



[see, e.g.. iWyrowski et al.ll2003l ). The frequencies in Table 1 represent the rest frequencies 
(for identifled lines) or laboratory frame frequencies, corrected using Visr = —23 km s"^ (for 
unidentifled lines). "U" in Table 1 indicates unidentifled lines. The radial velocities of the 
identifled lines in Table 1 are velocities at the intensity peak. No recombination lines lye 
within the observed frequency ranges. We inspected intensity maps of all detected molecular 
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lines, but we did not find any extended features. 

4. Discussion 

In tliis section, we focus on the properties of the ^^CO J = 6—5 hne. In particular, we 
discuss the structure probed by the ^^CO J = 6—5 hne, and the large missing flux of the 
low velocity component of the ^^CO J = 6—5 hne. 

A notable finding in the present observation is that the cavity structure, which has been 
found by the Sanchez Contreras's observation in the ^^CO J = 2—1 hne, is not clear in 
the ^^CO J = 6—5 maps, although the beam sizes are almost the same with the Sanchez 
Contreras's observation. An immediate interpretation for the missing cavity is that the 
emission of the ^^CO J = 6—5 line corresponding to the cavity wall is resolved out by the 
intcrfcrometry. However, we have to be careful in this interpretation, because the cavity wall 
probed by the ^^CO J — 2—1 line has exhibited thin structure (thickness~ 1"). 

In a careful second look of Figure 3, we reahze that features seen in velocity-channel 
maps of the ^^CO J = 6—5 seem to trace the south wall of the cavity even though it is 
not very clear. In Figure 3, for example, the map at —26 km shows an extension to 
the east-southeast direction from the nebula center, and the map at —16 km s~^ shows an 
extension to the west-southwest direction from the nebula center. In addition, in the lower 
panel of Figure 4, the intensity peak of the CO emission seems to be shifted (~ 0.8" to south) 
from the position of the continuum sources. These phenomenon, presumably, suggest that 
we detected only the south wall of the cavity in the ^^CO J = 6—5 hne, and that we did not 
detect the central dense core and the north wall of the cavity. 

However, we still have a couple of questions: (1) what is the source of the missing flux of 
the ^^CO J = 6—5 hne, and (2) why is only the north wall detectable in the ^^CO J = 6—5 
line? In terms of the expanding velocity and the angular size of the structure, only the 
spherical halo and the slow axial component could be a possible source of the missing flux. 
The emission from the spherical halo is expected to be partially resolved out if it is detectable 
in the ^^CO J = 6—5 line, because the spherical halo is largely extended (at least, in the 
^^CO J = 2—1 line). However, an important point is that we cannot expect the spherical 
halo to be highly excited enough to emit the ^^CO J = 6—5 line. The CO lines emitted from 
a circumstellar spherical halo are usually radiatively excited by infrared emission from the 
central star and dust component in the envelope (the energy input to the dust component 
originally comes from the central star. Therefore, we cannot expect a strong IR emission in 
the outer part of the envelope even if we assume the IR emission from the dust). Therefore, 
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only the cavity wall (slow axial component) may be the possible source of the missing flux, 
at least, among the known components in the CRL 618 molecular envelope. 



In the current scenario (iLee fc Sahail l2003l : ISanchez Contreras et al.l |200J) , the cavity 



wall is formed by an interaction between ambient material and the energetic high velocity 
jet. The shock caused by the energetic jet would be a possible reason for a high tempera- 
ture/density which would excite the ^^CO J = 6—5 line. In fact, the Sanchez Contreras's 
model, which has reasonably reproduced their observations in the CO and HC3N lines, has 
applied 200 K as the temperature of the cavity wall. This temperature is much higher than 
the excitation temperature of the ^^CO J = 6—5 line (~ lOOK). 

In our opinion, the following two reasons cause the absence of the north wall of the cavity. 
First, the temperature and/or density of the north wall is smaller than those of the south wall. 
Consequently, the north wall exhibits a weaker intensity in the ^^C O J = 6—5 line. In fact 



the n orth wall shows a weaker intensity in the ^^CO J = 2— 1 line (ISanchez Contreras et al. 



2004 ). Second, the different -uf-coverage between Sanchez Contreras's and the present ob- 
servations possibly causes the missing northern cavity wall. The rate of the missing flux 
strongly depends on the pattern of the wf-coverage. In some cases an elongated structure 
could be resolved out even if it shows thin structure. Presumably, the missing northern 
cavity is a combination of these two effects. The difference in the temperature/density be- 
tween the north and south walls, if it is real, leads to the fact that the physical conditions 
in the CRL 618 molecular envelope is not exactly axial symmetric. The southern part of 
the envelope (especially, in the innermost part) may exhibit a relatively higher tempera- 
ture/density as compared with the northern part. This non-axial symmetric distribution of 
the temperature/density might be caused by the central star (for example, binarity and/or 
time- variation of post-AGB wind, which are very poorly understood). 



5. Summary 

This paper has reported the results of an SMA observation of the proto-planetary nebula 
CRL 618 in the ^^CO J = 6—5 line. We also observed simultaneously in the ^^CO J = 2—1 
and ^^CO J = 2— 1 lines. In addition to the CO lines, we have detected a number of emission 
lines at 230 GHz. The main results are as follows: 

1. Flux of the high velocity component of the ^^CO J = 6—5 line is almost fully recovered, 
while we lost roughly 80% of the low velocity component flux. This suggests that the 
emission region of the low velocity component of the ^^CO J = 6—5 line is largely 
extended. 



- 11 - 



2. The existence of the cavity structure, which has been found in a previous observation 
in the ^^CO J = 2—1 hne, is not clear in the ^^CO J = 6—5 maps. Only the south 
wall of the cavity seems to be detected in the ^^CO J = 6—5 line. This result suggests 
that the physical condition of the molecular envelope of CRL 618 is not exactly axial 
symmetric. 

3. Continuum emission is detected both at 230 and 690 GHz. The flux of the 690 GHz 
continuum emission seems to be partially resolved out by the intcrfcrometry. This 
possibly suggests that dust emission significantly contaminates the 690 GHz continuum 
flux. 



We are grateful to Rob Christensen, Ahson Peck and the SMA team for making the 
690 GHz campaign possible. This research has been supported by the Academia Sinica 
Institute of Astronomy & Astrophysics and the Smithsonian Institute, and has made use of 
the SIMBAD and ADS databases. The authors thank Jinhua He, Holger Miieller and Shuro 
Takano for their help in the identification of the molecular lines. JN thanks Paul Ho for his 
constant encouragement. 
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Fig. 1. — S chematic view of the mol e cular envelope of CRL 618 made using examples from 
Figure 5 in ISanchez Contreras et al.l (120041 ) . The small arrows indicate the velocity field of 
the envelope. The numbers in parentheses correspond to the numbers of the components 
referred to in the text (see Section 1). 
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Fig. 2.— Spatially integrated spectra of the ^^CO J = 2-1, ^^CO J = 6-5 and ^^CO J = 2-1 
lines. The integrated area is a circle with an angular diameter of 15". For the ^^CO J = 6-5 
line, the spectrum is smoothed over every 5 km s~^. The other detected lines seen in the 
spectra are indicated by the name of the molecular species or "U" meaning unidentified line. 
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Fig. 3. — Velocity channel maps of the CO J = 6—5 line. The velocity channels are averaged 
over 10 km s~^ intervals. The contours start from a 3 o" level, and the levels are spaced every 
5 (7. The 1 a level corresponds to 4.42 x 10~^ Jy beam~^. The dashed contour corresponds 
to —3 a. The peak intensity corresponds to 28 a. The synthesized beam is indicated in 
the lower right panel. The background gray scale represents the intensity of the 690 GHz 
continuum emission. 
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Fig. 4. — Velocity integrated intensity maps of the CO J = 6—5 l ine superimposed on 
the HST WFPC2 H„+continuum image (ITrammell fc Goodrich! |2002| ). The red- and blue- 
contours represent the red- and blue-shifted parts of the high and low velocity components. 
The velocity ranges of the integration are given in the panels. The lowest contour corresponds 
to a 3 0" level, and the levels are spaced every 4 a; a 1 a level for the upper and lower panels 
are 1.2 x 10~^ Jy beam~^ and 6.9 x 10~^ Jy beam~^, respectively. The synthesized beam 
size is given in the lower-right corners. The dotted arrows, A and B, indicate the directions 
along which the p—v cuts shown in Figure 8 were taken. 
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Fig. 5. — Velocity channel maps of the ^"^CO J = 2—1 line. The velocity channels are 
averaged over 3 km s~^ intervals. The contours start from a 3 cr level, and the levels are 
spaced every 10 a. The 1 a level corresponds to 9.5 x 10~^ Jy beam""^. The dashed contour 
corresponds to —3 a. The intensity peak corresponds to 40 a. The background gray scale 
represents the intensity of the 1 mm continuum emission. The synthesized beam is indicated 
in the lower right panel; the outer unfilled and inner filled ellipses represent the beam sizes 
for the line and continuum maps, respectively. 
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Fig. 6. — Velocity channel maps of the low velocity component of the CO J = 2—1 line. The 
channel velocities are given in the upper left corners of each panel. The velocity channels 
were averaged over 1 km s~^ intervals. The contours start from a 3 a level, and the levels 
are spaced every 20 a. The dashed contour corresponds to —3 cr. The 1 a level corresponds 
to 1.03 X 10~^ Jy beam~^. The peak intensity corresponds to 120 a. The background gray 
scale represents the intensity of the 1 mm continuum emission. The synthesized beam is 
indicated in the lower right panel; the outer unfilled and inner filled ellipses represent the 
beam sizes for the line and continuum maps, respectively. 
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Fig. 7. — Velocity channel maps of the high velocity component of the CO J = 2—1 line. 
The velocity channels are averaged over 1 km s~^ intervals. The contours start from a. 3 a 
level, and the levels are spaced every 10 a. The dashed contour corresponds to —3 a. The 
1 a level, background gray scale and synthesized beam are the same as in Figure 6. 
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Fig. 8. — Position-velocity diagrams of the CO J = 6-5, CO J = 2-1 and ^^CO 
J = 2-1 lines superimposed o n the similar diagram of the ^^CO J = 2-1 line taken from 



Sanchez Contreras et al.l (120041 ) . "A" and "B" mean the direction of cut indicated in Figure 



4. The contours start at a 3 a level. A linear scale is used for the top panel, and the in- 
crement of the contours is every 1.75 Jy beam~^. A logarithmic scale is used for the middle 
and bottom panels. The dashed horizontal lines represent the origin of the offset axes. 
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Fig. 9. — Spatially integrated spectra in the 1 mm band. The integrated area is the same 
as Figure 1. The horizontal axis represents the laboratory frame frequency, corrected using 
Visr — —23 km s~^. The vertical sohd bars indicate the rest frequencies of identified lines. 
"U" means unidentified hne. 
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Table 1. Line Detections 



Frequency 
(MHz) 


Species 




Transition 




isr 

(km s~^) 


s 

(jy km s 


219675.114 


HC3N 


J = 


24-23 i/7 = 2 Z = 





-22.0 


26.3 


219707.349 


HC3N 


J = 


24-23 1/7 = 2 Z = 


2e 


-22.1 


25.4 


219741.866 


HC3N 


J = 


24-23 i/7 = 2 Z = 


2/ 


-23.1 


25.1 


219866.826 


U 








— 


1.37 


219892.538 


U 








— 


1.13 


219934.820 


C^^N 


N ^ 


:2-l J = 5/2-3/2 




-22.4 


4.89 


220050.478 


U 








— 


2.89 


220070.185 


HC3N 


J = 


24-23 1/7 = 3/ = 


le 


-20.2 


10.1 


220181.972 


u 








— 


4.87 


220278.206 


u 








— 


2.54 


220332.566 


u 








— 


2.54 


220398.681 




J = 


2-1 




-10.0 


310.0 


220699.688 


HC3N 


J = 


24-23 z/7 = 3 Z = 


1/ 


-21.2 


17.4 


220709.024 


CH3CN 


Jk-- 


= 123-1 I3 




-22.5 


14.9 


220730.266 


CH3CN 


Jk-- 


= I22-II2 




-21.6 


16.5 


220743.015 


CH3CN 


Jk-- 


= 12i-lli 




-22.3 


— 


220747.265 


CH3CN 


Jk-- 


= 12o-llo 




-25.5 


— 


220868.094 


u 








— 


2.60 


220932.401 


HC5N 


J = 


83-82 1/ = 




-23.5 


15.8 


220957.716 


u 











3.28 


220979.090 


/-HC4N 


J = 


48-47 




-22.9 


3.08 


221020.158 


u 










3.34 


221090.680 


u 










8.77 


221246.240 


Hi^CCCN 


J = 


25-24 z/7 = 1 Z = 


le 


-18.3 


19.8 


221295.423 


HC5N 


J = 


83-82 = 1 / = 


= 1/ 


-23.3 


5.15 


221486.033 


HC5N 


J = 


83-82 un^ll^ 


= le 


-22.8 


2.42 


229589.088 


HC3N 


J = 


25-24 i/7 = 3 Z = 


3/ 


-19.9 


16.5 


229772.876 


H^^CCCN 


J = 


26-25 1/7 = 1 Z = 


le 


-23.7 


6.61 


229891.277 


HC3N 


J = 


25-24 1/7 = 3 Z = 


le 


-23.6 


9.46 


230074.654 


U 










2.19 
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Table 1 — Continued 



Frequency 


Species 




Transition 




S 


(MHz) 








(km s-i) 


(jy km s-1 ) 


230093.224 


Hi^CCCN 


J = 


26-25 ur^ll^lf 


-18.8 


14.6 


230180.393 


/-HC4N 


J = 


50-49 


-17.4 


4.04 


230299.796 


U 








5.81 


230538.000 


CO 


J = 


2-1 


-20.0 


1541.4 


690552.098 




J = 


8-7 


-19.6 


346.4 


691473.076 


CO 


J = 


6-5 


-22.0 


4647.8 



